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Abstract 

Recent theories have proposed a variety of massive particles, like the moduli, whose abundance or 
decay endangers standard cosmological results. To dilute them, thermal inflation has been proposed, 
with its own massive scalar flaton held which goes on decaying in the era of Mev-scale temperatures. 
In this paper, the effect of late-time entropy production on neutrino decoupling during such reheating 
is investigated by including a term, arising from the rate of entropy production due to scalar decay, in 
the Boltzmann equation for the neutrino number density. The effect on the decoupling temperature 
of massless neutrinos is studied. It is found that a lower bound to the scalar decay rate constant 
can be set at 10~ 22 Gev. 
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1 Introduction 



Various massive fields like the gravitino, the Polonyi, the moduli, and the dilaton figure in 
supersymmetric and string theory models [fij. The corresponding particles are long-lived. 
These fields pose quite serious cosmological problems. If they decay during baryogenesis or 
nucleosynthesis(BBN), the baryon to photon ratio may be greatly diluted in the first case 
and nuclear abundances distorted in the second. If they are stable, there is a problem of 
over-abundance. People have studied the problem in detail and proposed various ranges of 
masses and decay rates to minimise damage to standard cosmology || ||, [|. 

An interesting proposal is that of thermal inflation [||, H. Quite apart from primordial 
inflation |7|], a scalar field, called the flaton, is used to generate inflation at late times so 
that the potentially dangerous fields mentioned above may be diluted away. Typically ||, 
the inflation starts as the thermal energy density falls below the flaton potential at the 
origin at a temperature of about 10 7 Gev, and stops when the flaton vev reaches 10 12 — 10 14 
Gev at a temperature ~ the flaton mass, which may be taken to be < 10 3 Gev. Such a 
flaton field will go on decaying into the Mev-scale temperature era, and the parameters in 
thermal inflation must be such that the entropy generated due to flaton decay may not 
affect nucleosynthesis. 

Recently, the effect of the decay of a massive particle, like the flaton, on nucleosynthe- 



sis has been subjected to detailed investigation [||, [l0|, [Tl[]. As the parameters and decay 



schemes of such particles are yet to be generally agreed upon, it is useful to address the 



general problem of the effect of late-time entropy production, during reheating |T^, 13, 14], 
in the context of BBN. In the present work, neutrino decoupling is studied in the presence 
of such late-time entropy production, neutrino decoupling temperature being one of the 
key parameters determining BB nucleosynthesis. Bernstein's method of "pseudo chemical 
potentials" [15] is used to introduce the entropy production rate directly into the Boltzmann 



equation for the neutrino number density. The decoupling temperature is estimated in the 
presence of this entropy generation. 

The paper is organised as follows. Section 1 is this Introduction. In Section 2, the 
late-time development of a universe, once dominated by a scalar field, is set out with 
some modification of the usual formalism. In Section 3, the method of "pseudo chemical 
potentials" is applied to modify the Boltzmann equation for the neutrino number density 
and incorporate into it the late-time entropy production rate directly. In Section 4, the 
effect of late-time entropy production on neutrino decoupling is studied, using the results 
of the previous sections. In Section 5, conclusions are stated. 

2 Reheating at Late Times 

In this section, a two-component universe is considered, made up of the scalar field energy 



density p<p and a radiation energy density pr. It is assumed [16, 0], further, that either the 



scalar field parameters have not supported resonance, or, in any case, the times considered 



are late enough to be beyond a possible preheating [17] interlude. 
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Taking the scalar decay rate constant to be T 13], the equation for the evolution of the 
scalar energy density is 

d 

— ptj, + 3H P(j , = -Tp,p. (1) 
Defining = a 3 and R = a 4 /^, where a is the scale factor, (||) becomes 

$ = -r$. (2) 



From 



one obtains 



d o d 3 

^[a (p<* + Pfl)J +PR-Q t a =0 ' 



i? = ar$. (3) 



The Friedmann equation is 



5. 



/f 2 = -^a + fl), with ^ = ^p-. (4) 



2.1 Incomplete ^-domination 

In this subsection, it is assumed that R <C $a, so that the square and higher powers of 
Rj ($a) can still be neglected at these times. (|j) gives 

H = --H 2 (l H — - — — -). 

2 v 3($a + R) J 

With the R <C &a asssumption, one can write 

H^--H 2 {1 + — ). (5) 
2 v 3$a v ; 

For full matter domination, at yet earlier times, R is neglected compared to <3?a, and 

2 

leading to t = ^j. In the era of incomplete ^-domination, the correction term on the RHS 
of (H) is not neglected, and may be evaluated to the approximation 

R-R I = ^(a-a I ), (6) 
da 

where ti refers to some initial epoch such that a » a;. Also, if it is supposed that the scalar 
decay produces sufficiently copious radiation, Ri <C R. As a correction term is being dealt 
with, these approximations should not cause much deviation from the actual evolution. 
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Then, for t sufficiently later than tj, but within the regime under consideration, one may 
write, in the correction term on the RHS of (@), 

R « ^a, (7) 

using (H) and (g). (|B|) becomes 



Now, evolution will be described by the new variable 

r 

H 



r 

Z = T7, (9) 



instead of the time t. (01) can be written as 



1 + x), (11) 



a 3 



using ([?]) in the correction term on the RHS of (jlOj) . 
Then, (g) leads to the evolution equation for $ 

(12) 

ax 3 + x 

with the solution 

$ = ^ V. , i d » (13) 



in the approximation i « 1. To obtain R, (|l^) is substituted in @, a being calculated 
from (11), using (|). The result is 

dR 2.K.2 2 4 (1 + 1^)1(1 + s )| 
dx 3 T (i+i x )t- 

Integration will give R. 

To compare with earlier work, at early times, one goes to the approximation x,xj -C 1. 
( |l4l ) becomes 

On integration, taking i?/ <C -R, 



2, #,2,1, 5 I, 

5 



2 r 



i 



<3?f (a 2 - af ), 



which agrees with [16]. 
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2.2 Incomplete Radiation domination 

Now, one may consider even later times when $a <C R, such that <&a/R cannot be neglected, 
but its higher powers can. The Friedmann equation reads 



K 2 „, $a 

U" = 

which leads to 



H 2 = ^ T R(1 + ^), (15) 
a it 



H = -2H 2 [l- — — - ]. (16) 
L 4(<I>a + i?) J v ; 



If, well into this epoch, the correction term on the RHS of (|16|) is neglected, the full radiation 
domination relations are found: 

H = — , and 

a = At%, (17) 

A being a constant. 

(Eh has, as solution, a falling exponential in i, viz. $ ~ e~ r *. The new evolution variable 
x = T/H is now sought to be introduced in place of t. To do this, instead of taking the 
falling exponential in t directly, suitable approximations to the correction terms on the RHS 
of fllq) and fll6j), are first worked out. Let to be a sufficiently late epoch, when $ = $0 ~ 0. 
Then, for use only in the correction terms, one takes 

$ _ $0 = $ (_)_ $ (_) = _ k( ___). 
Neglecting l/*o compared to $, respectively, an approximation 

**f, (18) 

will be used only in the correction terms, i.e. in the correction terms, the falling exponential 
will be approximated by a rectangular hyperbola. B is a constant. Using this approximation 
in the correction terms, the final form of will be found below to be a falling exponential 
in (l/2)sc, with a pre-exponential correction. 

Next, a similar approximation is considered for R. It ought to be mentioned that R 
refers to the total radiation present, and not only to that produced by decay. However, the 
change in R is due to <j) decay and consequent entropy production . In the absence of this 
decay, R = 0. 

Using ( |l7l) and (|lg| ) in (j3|), and, integrating, one obtains, for use only in the correction 
terms, 

R-R E ^ 2ABT(t? - if). 

If tE is sufficiently early compared to t, though within the regime under consideration, and 
there is sufficiently copious radiation production since tE, it is sufficient to take 

R « 2ABTt% 
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in the correction terms. ( |i~7| ) and (|T^) are now used to give, in the correction terms, 

— - — , once again, as m 7), 

<3?a H 

= x. (19) 
The approximation is now x S> 1. (|l5| ) and ( |l6|) become 



and, # = -2# 



K 2 1 

"'1 + 

x + 



H 2 = ^R(l + -), (20) 
a* x 

3 

2^ 1 4 



X + l 

This last equation and (||) give 

$ = $ E (l^±^)| e -|(— a) (21) 
Ax + 3 

One can see that this solution is dominated by the falling exponential. The pre-exponential 
factor arises because radiation domination is not yet complete, and (l/2)x is not quite Tt. 
For sufficiently late times t,tE, one has x,xe ^> 1, (20) gives H = l/(2t), in the usual way, 
and reduces to 



which agrees with [12]. 
From @, 

a»(^)3fli(l + i)i. (22) 
Now, on integrating (||), using (|2~l|) and (|22]), one has 



i?4 = 



4 

oo 



3 »ii . i . 3 N i i_ f 00 3 i_, . 

-^Y ) ^ { - Xe + a } E I x ~ 8e ~" ' (23) 

for i > 1, where i?oo is the value of R after decay is finished. The integral in (^) is an 
incomplete Gamma function, and can be numerically evaluated for any value of x. 
The temperature is defined from 

R = g *—a A T\ (24) 
y 30 v ' 

where g* is the effective number of relativistic degrees of freedom. 
The entropy production originates from [18] 

dS^ = (25) 
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3 Entropy Production and the Boltzmann Equation for the 
Number Density with Scalar Decay 

The effect of entropy production on decoupling is studied here with reference to the par- 
ticular problem of the electron neutrino, u, decoupling from a thermal bath of electrons, 
positrons, and photons, in the presence of a decaying scalar field. It is assumed that the 
scalar field has a very small branching ratio into neutrinos |9]]. So, the main process which 
changes the number density n of the v neutrinos is 

v(k)+P(k)—>F(p)+F(p), (26) 

where J-,F are fermions, and, k is the energy-momentum 4- vector (Ek,k). Apart from 
e~,e + , the relevant fermions, at the epoch of v decoupling, are the neutrinos of the other 
families. 

If there is no other process like scalar decay, it is usual to assume [[0| that, in the 
absence of Fermi degeneracy, the distribution functions may be written as 

f u (k) = f(k) = e-^W"^* 

f e (p) = 

f P (k) = f'(k) = e -°'(*)-/^s 

fe+(P) = e-^, (27) 

where T = 4 is the temperature of the thermal bath, and a' is a time-dependent "pseudo 
chemical potential" , introduced to take into account the departure of the decoupling neu- 
trinos from equilibrium. (In (27), it has been assumed that T is the electron.) 
Decoupling is, in this case, governed by the integrated Boltzmann equation 

h + 3Hn = - < a\v\ > (r? - n 2 EQ ), (28) 

where ueq is the equilibrium number density, 

nEQ = I J^ e '" Ek = r3/7r2 < (29) 

n^n EQ e- a '[15], (30) 

and 

< a \v\ >= ^— [ dRdKe-^e-^I, (31) 

n EQ J 

I being the invariant integral 

/ = J dPdP(2ir) 4 5 4 (p + p-k-k)\M\ 2 . (32) 

g is the number of spin degrees of freedom, dK = gd?k/[(2TtY2Ek\ etc., and \Ai\ 2 the 
spin-averaged square of the modulus of the relevant matrix element. 
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Considering the process (p6|), and, assuming CP-invariance and the absence of Fermi 



degeneracy, the Boltzmann equation for the neutrino number density n may be written[15 



^§- t (a 3 n) = -]T J dKdKdPdP(2^)H\p + p-k-k)\M T \ 2 (f'f' - //"). (33) 



/' corresponds to the v neutrino and / to the fermion T . 
Then, 

-^(a 3 n) = - J ' dKdKdPdP{2ir) 4 5\p + p-k-k)\M\ 2 (fT - //) 

- Y,( dKdKdPdP(2TrfS\p + p-k-k)\Mi\\f'f'-fifi), (34) 

i 

where, on the RHS, the first term refers to the process 

v + v — >e~ + e + , (35) 

and the second to the processes 

v + v — ► Vi + i = fx, t, (36) 

fi being the U{ distribution function. 

The injection of entropy will cause changes in the distribution functions. One way to 
tackle the situation is to use the Boltzmann equation for the v distribution function, instead 
of the integrated Bolzmann equation for the number density, include u+e~{e + ) — > z/+e~(e + ) 
processes, and proceed numerically ||. However, the integrated Boltzmann equation for 
the number density will be used here to study decoupling at late times, when the entropy 
generation rate is small, in the following way. 

It has been assumed^], in accordance with current models, that the <fi does not decay into 
neutrinos. But, ref.|| finds that, for reheat temperatures < 7Mev, the neutrino distribution 
function is materially affected by the decay. How can entropy injected into the e, e + , 7 sector 
affect neutrino distribution function and, hence, neutrino density? The transfer of entropy 
must occur via process (|35|) , and this must lead to extra terms on the RHS of ( |28|) if n is to 
be affected. From (|33[), this can only occur if /, the electron distribution function, changes 
due to entropy injection. 

The electromagnetic interactions, faster than the expansion, will tend quickly to ther- 
malise the e,e + ,7 (i.e. establish kinetic and chemical equilibrium in this sector). But 



recent workjlE, has shown that this thermalisation of <\> decay products is not nearly 
instantaneous. A massive (ft, as it goes on decaying, continually injects highly energetic light 
particles, with energies much above the thermal mean. There is no possibility of inverse 
decay and so the term from scalar decay in the Botzmann equation for, say, the electron 
distribution function opposes the equilibriating effect of the electromagnetic interaction. As 
long as the decay does not become negligible, complete thermalisation of the distribution 
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functions of the decay products is not assured. How quickly complete thermalisation will 
occur, i.e. the epoch of thermalisation, depends, not only on the relative rapidity of the 
electromagnetic interactions and the expansion, but also on the parameters of cf). In fact, the 
requirement of thermalisation before a specific epoch (e.g. BBN) has been used to derive 
bounds on these parameters ]i~9|, p0|, 11]. 



The thermalisation of rapidly interacting, light decay products, of a heavy boson was 
considered in ref.|l5|]. There, the progress, in time, of the light particle distribution function 
g was directly considered. A solution g = e~ a ~^ E of the Boltzmann equation was mooted, 
where a was a time-dependent "pseudo chemical potential", although, the actual calculation 
was done with a 5 function distribution, because the problem was simplified by assuming 
that all the heavy bosons decayed at one instant. Following this lead, it will be supposed, in 
this paper, that the injection of entropy into the bath, due to scalar decay, can be taken care 
of by introducing a small time-dependent "pseudo-chemical potential" a(t) into the electron 
and positron distribution functions, in addition to the potential a'(t) already introduced in 
the v distribution functions. 

Of course, the electromagnetic interactions, being so fast, will cause a to be very small. 



In effect, then, the scheme in (27) will be changed, in the absence of Fermi degeneracy, 



to 



fu(k) = f'(k) 
fe(p) = f(p) = e~ a -^ 

f = e- a - pE f. (37) 

The entropy density Sj, contributed by a particle j with a distribution function hj(q), 
may be defined, in the absence of degeneracy, as 



Following [15], the covariant divergence of the entropy density current for a component j, 



in a process like that described by (|35|), can be shown to be 

1 d , 3^ f d 3 q Cj. , 

^di {Sja) ™-J * (2^^' 



where Cj is the relevant collision integral as defined b, f\in [18]. For example, for the electron, 

i__a 

a 3 dt 



' d -(s e a 3 ) = - J dKdKdPdP(2n) 4 5 4 (p + p-k-k)\M\ 2 



x Inf(ff'-ff). 



3 This differs from the definition in ]f5(] by a factor of Ej 
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It can immediately be seen that the covariant divergence of the total entropy density current 
for the process (^5|) is, assuming energy conservation, 



X^ristota 3 ) = [ dKdKdPdP(27r) i 5 i (p + p-k-k)\M\ 2 
a 6 ot J 

x2(a'-a)(e- 2tt -e- 2a ')e-^^) 

J dKdKdPdP(2ir) 4 5 4 (p + p-k-k)\M\ 2 

x2a'(l-e- 2a ')e- /3(£fc+%) 

+ J dKdKdPdP(2ir) A 5 A (p + p-k-k)\M\ 2 

x (-2a)(l - e~ 2a ' + 2a')e-^ Ek+E ^ . (38) 

Here, a, being the effect of entropy injection, is assumed to be small at late times. Squares 
and higher powers of a have been neglected. 

The first term on the RHS of ( |38|) is the covariant divergence of the entropy density 
current for a = 0, and the second term is proportional to a. So, it can be concluded that 
the second term measures, approximately, to first order in a, the contribution of eft— decay 
to the covariant divergence of the entropy density current, corresponding to process (|35|), 
viz. 

4^ = -2a(l - e~ 2a ' + 2a') / dKdKdPdP{2Tif5 i {p + p-k-k)\M\ 2 
a 6 ot J 

xe -P(E k +E- k ) ^ (39) 

where is defined in (|25| ) and e is the fraction of the entropy current from 0-decay which 
contributes to process (p5[). 

Again, in the presence of <f>— decay, 

i-^-(a 3 n) = - ( dKdKdPdP{2ir) A 5 A (p + p-k-k)\M\ 2 
a 6 ot J 

X ( e -2a> _ e -2* )e -P(E k+ E k ) 

- E/ dKdKdPdP(27T) 4 5 4 (p + p-k - k)\Mi\ 2 (ff' - fifi) (40) 

i 

= -(n 2 - n 2 EQ ) < a\v\ > 

- 2a J dKdKdPdP{2it) A 5 A {p + p-k- k)\M\ 2 e~^ Ek+E ^ 

- Yl[ dK dKdPdP{2^) A S\p + p-k-k)\Mi\ 2 (f'f'-fifi) (41) 

i 

to first order in a, using (|3C|), (|3l|), ( ^7| ) and (p4]). < > refers to (35). 

Now, the finer point that, in this epoch, the electron neutrino has both charged and 
neutral current interactions, while the other neutrinos have only neutral current interactions 
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(muons and taons having already decoupled), will be neglected, and it will be assumed that 
the charged and neutral current interactions fall out of equilibrium together. Also, neutrinos 
of all types will be assumed to be massless. Then, apart from their contributions to pr as 
separate species, the neutrinos of the different types may be treated as identical, and it may 
be supposed that they decouple together. In this approximation, the difference between /' 
and fi may be neglected, and the terms within the summation sign dropped, in comparison 
to the other terms on the RHS of 

So, to first order in a, the integrated Boltzmann equation can be written, in the presence 
of (ft— decay, as 

1 d o 

— ? — (a n) = h + 3Hn 
a A at 

= -(n 2 -nl Q )<a\v\> + {1 _J al + 2a , ) ^, (42) 
from (|39|) and (|4l|). ([42]) will replace (|28|) in the presence of scalar decay. 



4 Effect of Scalar Decay on Neutrino Decoupling 



Decoupling will be supposed to set in when the RHS of (|42| ) becomes less than 3Hn. To 
analyse the two terms on the RHS, radiation domination is first assumed. A change of 
variables is introduced: 



I 3 3 

— ,Y = nx2,Y EQ = n EQ x2. 



Also, ( pop gives, in the era of neutrino decoupling, 



TT 4.461 x 10~ 19 m2 / Is 1 
H= T 2 1 + - 2. 

Gev x 



(43) 



(44) 



taking g* = 10.75. The LHS of @, then, becomes (2r/xa )dY/dx. Near the start of the 
decoupling regime, one can put 

Y = Y EQ + A, 

where A is very small, and it is possible [18, 22] to set 



— — = 0, and, 
dx 

dY dY EQ 
dx dx 



(45) 



(g|), @ and (|9|) are used to find that dY EQ /dx > 0, i.e. the LHS of (f42|) is positive. Use 
of (pp|) shows that each term on the RHS of (|42|) is of the same sign. So, for the equation 
to hold, 

a' > 0, (46) 

and each of the terms on the RHS of (^) is positive. This means, that when decoupling 
sets in, each term must be separately less than 3Hn. 
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TV, f A/TcoA 

i £) \ivi ev ) 




0.75 


4.33275 


1 


3.4705 


1.25 


2.8038 


1.5 


2.264 


1.75 


1.8174 


2 


1.4470 


2.25 


1.1438 


2.5 


0.901 


2.75 


0.7109 


3 


0.5644 



Table I: Neutrino a' values at decoupling for different decoupling temperatures 

without scalar decay 



The value of a' chosen will determine the temperature of decoupling. This choice is 
fixed as follows. In the absence of (ft— decay, there is only the first term on the RHS of (|42[). 
By setting this term equal to 3Hn, a set of values for a' is obtained for neutrino decoupling 
temperatures To between 0.75 and 3 Mev, a range which more or less safeguards BBN. It 
may be mentioned that, at present, observed 4 He abundance values, Y p , fall between 0.234 



1 23(1 and 0.244 | 24j . A rough estimate, as indicated in [18], gives corresponding decoupling 
temperatures ~ 1 Mev. 

Next, for these values of a', the full RHS of (|42| ) is put equal to 3Hn, and the decoupling 
temperatures, corresponding to different values of the scalar decay constant T, are worked 
out. 

The first step uses 

(n 2 - nL) < a\v\ > 

- ~ ikr^- £ L (47) 

As already discussed, in the evaluation of < a\v\ >, it is enough to consider the process 
(|35|), in the s-channel with Z exchange and the t-channel with W exchange. Evaluating 
(|3|), and, then, ©, using (g|), 

<aM> = -G 2 F [{C Ve + lf + {C A e + lf]T 2 

7T 

4.112 x lO" 1 " n ,„ . 

" a* T - (48) 

The mass of the electron has been neglected. Gf is the Fermi constant. 

Taking H = 4,461 g e ^° — T 2 , in the absence of </>-decay, and using (p9[), (p0[), and (P 
one obtains the values of a' shown in Table I. 



Now, the second term on the RHS of (42), i.e. the contribution of scalar decay, is 
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considered. Using (f2l|), this term becomes 



e 1 dSj, e $_ 

(1 - e~ 2a ' + 2a')a?~df ~ ~ (1 - e~ 2a ' + 2a') Ta 3 



4ic+3 



'1 - e" 2 "' + 2a') Ta 3 



(49) 



This term, for x >> 1, will be dominated by the exponential. So, rather drastic assumptions 
will be made to estimate the pre-exponential, in the absence of phenomenological informa- 
tion about xe and &e- xe refers to any epoch sufficiently early in the regime of incomplete 
radiation domination, and is set equal to 1, so that &e = Re/o-e- This is quite an approxi- 
mation, because (21) is faithful to ~ e~ Tt provided xe,x 3> 1. With this approximation, 
Te is, then, estimated by setting x = 1 in ([44|). Next, one puts R E = (tt 2 / 30) g* E a%T E and 
9*SE a E^E = ds^T 3 . This last is the weakest assumption because it is only true if there 
is conservation of entropy, with a oc T~ l behaviour. But, entropy is being produced, and, 



in fact, in the regime of $ domination, T a [12]. In the early part of the regime of 
incomplete radiation domination, too, one expects a deviation from a oc T _1 . However, the 
expectation is that the exponential dominates the entropy current and its effect on v decou- 
pling, and one can take the view that the pre-exponential is being estimated only to order 
of magnitude. Neglect of the deviation from a oc T _1 behaviour in the pre-exponential, and 
the other approximations, should not make too much difference to an order of magnitude 
calculation of lower bounds on the scalar decay constant. 

With all these assumptions, the condition of neutrino decoupling in the presence of 
scalar decay becomes, using (gj), 

( n2 ~ n 2 EO ) < a \ v \ > . e IS* 



+ 



3Hn (1 - e~ 2a ' + 2a') a 3 3Hn 

208.945 sinha'vl el0.4927xt(l + x)*e~^ x ~ 1) 



x^(x + l)^ ' (x + |)f(l - e~ 2a ' + 2a')e 



where 



r 







W~ 22 Gev' 

It must be remembered that x > 1 in (|50|), corresponding to radiation domination. 



4.1 Numerical Results 

Three values of a' are used from Table I, namely, those corresponding to decoupling in the 

absence of scalar decay at Tjj = 1,2, 3Mev. In each case, the RHS of ( |50| ) is put equal 

to 1, and the resulting equation solved for x = xd, taking different values of Tq. The 

corresponding decoupling temperatures T' D are found from (|44|), using x = F/H to rewrite 

it as 

T = 14.9718^ 

Mev (X 2 +X)2 
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T D =lMev 




'l' D =2Mev 




Tn—SMev 




«'=3.4705 






a'=1.4470 




a'=0.5644 


r 


X D 


rpl 


r 


X D 


rpl 


r 


X D 


rpl 


100 


22414 


1.00 


100 


5603 


2.00 


100 


2490 


3.00 


10 


2241 


1.00 


10 


560 


2.00 


10 


248 


3.00 


1 


223 


1.00 


1 


55.2 


2.01 


1 


24.2 


3.01 


0.1 


22.6 


0.985 


0.25 


16.3 


1.83 


0.5 


15.4 


2.65 


0.02 


16.5 


0.513 


0.1 


13.7 


1.26 


0.1 


12.9 


1.30 








0.02 


13.0 


0.576 


0.02 


12.7 


0.583 



Table II: Decoupling Temperatures T' D in Mev for different values of 



the Scalar Decay parameter r = T/(10 22 Gev) 



e would be 1 if the entire entropy current of decay contributed to process fl35|). There 
are two points here. The decay modes of </> are unknown, and it is assumed that decay 
products, more massive than the electron, rapidly decouple, and almost all the entropy 
current of decay ends up first in the e,e + ,7 sector and, then, goes on to contribute to the 
entropy current of process ( |35| ) via a non-zero a. Also, a/fl means that there will be a 
small contribution of the entropy current from (ft— decay to the e, e + , 7 processes also. But, 
because a is so small, e has been put equal to 1. 

The results are shown in Table II. The results indicate that, for values of the scalar decay 
constant V > lCT 22 Get>, the neutrino decoupling temperature is not appreciably affected 
by scalar decay. 

It is necessary to check what happens if there is matter domination. Taking xj « 1 
in (13), extrapolating this equation to $ = at x = :c£; = l, and evaluating with the 



same set of assumptions as in the case of radiation domination, the criterion of decoupling 
will become, in place of 



(n 2 ~ n% Q ) < a\v\ > e 1 S 4 



3Hn (1 - e- 2 °' + 2a') a 3 3Hn 

208.945 sinha! rf 17.3934x 5 / 4 (l + x) 1 / 4 



3 

1 7 3Q3zLt. 5 / 4 M J- tW 4 

< 1. (52) 



a*(x + l)i (l-e- 2 «' + 2a')e-«'(l + ix) 

Matter domination implies that x < 1 in (|5^). Writing the two terms on the RHS of ( |52[ ) 
as A and B, A + B < 1. An examination of B will show that as the temporal variable x 
increases, B increases. So, if B predominates in an epoch, there can be no decoupling in 
that epoch. A little numerical work shows that for those values of T (greater than a critical 
value which depends on a') for which A predominates over B, A+B does not fall below 1, 
for values of x < 1 (matter domination ). So, no decoupling is possible in the epoch of 
matter domination. 

The numerical work shows, therefore, that for neutrino decoupling to proceed without 
significant change, in the presence of scalar decay, there must be radiation domination, in 
the sense that x = pr/ Pcf, > 1, and the scalar decay constant T must be > 10~ 22 Gev . 
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This corresponds to reheating temperatures > 8.5Mev , taking for the reheating tem- 
perature, the definition of ref.@ : 



T R . T ,1 2.4x 10 18 

0.554(- 



2 



Gev Gev y Gev 

In it was found that the neutrino distribution function was distorted, and the effective 
number of neutrino types, N e ff, started to decrease below three, as the reheating temper- 
ature fell below 7 Mev, and standard BBN was endangered. So, the present results show 
broad agreement with those of || on the question of the minimum reheating temperature 
which safeguards standard BBN. 



5 Conclusions 

The effect of entropy injection due to scalar decay on neutrino decoupling has been studied 
here by introducing the entropy production rate term directly into the Boltzmann equation 
for the neutrino number density. The method adopted was to introduce a small "pseudo 
chemical potential" a into the electron distribution function, in addition to the standard 
introduction of such a potential in the neutrino distribution function. As the electromag- 
netic interactions tend to thermalise the electron distribution function, this a is bound to 
be very small. Its function is to transmit the entropy current, arising from scalar decay, 
from the e, e + , 7 sector, to the neutrino sector through processes like e~ + e + — > v + v. 

The conclusion drawn regarding the condition of validity of standard BBN from this 
study of neutrino decoupling, viz. reheating temperature> 8.5 Mev, agrees broadly with 
the conclusions drawn in the literature from calculation of the form of the neutrino distri- 
bution function and the effective number of neutrino types, by numerical integration of the 
Boltzmann equation for the neutrino distribution function^. 

It is found that, in the presence of entropy injection due to scalar decay, the universe 
must not be matter dominated(p0 > ppt) in the neutrino decoupling epoch, as in that case 
there is no decoupling; rather, the universe must have attained radiation domination in the 
sense that p^ < pr. Consideration of the effect of scalar decay on decoupling temperatures 
in the regime of radiation domination leads to a lower bound of about 10 -22 Gev on the 
scalar decay constant V. 
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